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A theoretical investigation of the effects of elastic coherency on the thermodynamics, kinetics, and 
morphology of intercalation in single LiFeP04 nanoparticles yields new insights into this important 
battery material. Anisotropic elastic stiffness and misfit strains lead to the unexpected prediction 
that low-energy phase boundaries occur along {101} planes, while conflicting reports of phase bound- 
ary orientations are resolved by a partial loss of coherency in the {100} direction. Elastic relaxation 
near surfaces leads to the formation of a striped morphology, whose characteristic length scale is 
predicted by the model and yields an estimate of the interfacial energy. The effects of coherency 
strain on solubility and galvanostatic discharge are studied with a reaction-limited phase-field model, 
which quantitatively captures the influence of misfit strain, particle size, and temperature on solu- 
bility seen in experiments. Coherency strain strongly suppresses phase separation during discharge, 
which enhances rate capability and extends cycle life. The effects of elevated temperature and the 
feasibility of nucleation are considered in the context of multi-particle cathodes. 



Lithium iron phosphate (LiFeP04) emerged as 
an important high-rate cathode material for recharge- 
able batteries [1 and is unique because of its strongly 
anisotropic diffusivity [H [3] , its strong elastic anisotropy 
[T, and its tendency to phase-separate into lithium-rich 
and lithium-poor phases [5-9 . Despite a few conclusive 
observations of phase boundaries in chemically delithi- 
ated LiFeP04 nanoparticles [ 1QHT2] , the general consen- 
sus has been that phase boundaries always form during 
electrochemical discharge, thereby limiting battery per- 
formance Oils]. However, this limitation is inconsistent 
with dramatic rate improvements resulting from smaller 
nanoparticles, doping [14 , and surface coatings [15]. 

The feasibility of phase boundary formation has re- 
cently been challenged by both phase-field models [161 [17] 
and ab-initio calculations [18 . In a companion paper 
[17], we demonstrate that high discharge currents can 
suppress phase separation in reaction-limited nanoparti- 
cles, so that the spinodal is a dynamic property of in- 
tercalation systems. In this paper we consider the ad- 
ditional effect of elastic coherency strain and find that 
it leads to a quantitatively accurate phase-field descrip- 
tion of LixFeP04 that is useful for interpreting experi- 
mental data. Via mathematical analysis and numerical 
simulations of galvanostatic discharge, we conclude that 
coherency strain strongly suppresses phase separation, 
leading to better battery performance and improved me- 
chanical durability. 

Coherency strain arises when molar volume is a func- 
tion of composition, i.e. due to the difference in lattice 
parameters between FeP04 LiFeP04. Two-phase 
systems with identical crystal structure and small mis- 
fit strains generally form coherent interfaces [TO', "2^. 
It has been suggested that Li^FeP04 i"etains coherency 
throughout nucleation and growth [9], and in-situ obser- 
vation of crystalline material during battery operation 
supports this prediction [lU [22] . 

As we show, the observation of aligned phase bound- 
aries and striped morphologies in Li^FeP04 P^OW^] pro- 



vide conclusive evidence of coherency strain. Further- 
more, it is necessary to consider the fully anisotropic 
elastic constants and misfit strain to interpret experi- 
ments. Simplified elastic analysis has led to the con- 
clusion that {100} is always the preferred orientation 
[23^ 24 , although {101} phase boundaries are sometimes 
observed [TTJ [12] . Our fully anisotropic analysis predicts 
that {101} is the low-energy orientation, and we attribute 
the observation of {100} boundaries to a partial loss of 
coherency resulting from dislocations (or cracks). 

The origin of striped morphologies [TO^ has not been 
satisfactorily explained. It has been suggested that they 
result from the characteristic wavelength of spinodal de- 
composition [12 , although it is not clear why the insta- 
bility would be frozen in this state. We show instead 
that stripes represent the stable equilibrium state of fi- 
nite size particles, and predict that the spacing scales 
with the square root of particle size. As a result, we are 
able to extract the interfacial energy from experimental 
micrographs. 

The reported solubilities of Li in FeP04 v^iry signifi- 
cantly and depend on particle size and temperature [8]- 
|T0l|2T| 23 , 25 - 28 . These differences in solubility can now 
be explained in light of coherency strain. Phase field cal- 
culations of solubility as a function of particle size and 
temperature are able to fit experimental data with just 
two parameters. 

An alternative to coherent phase separation is for en- 
tire particles to remain homogeneous and form a mo- 
saic pattern, with some particles existing at low con- 
centration and others at high concentration [21] |29l [30] . 
This scenario is energetically favorable since there is no 
phase boundary energy or change in solubility, but it re- 
quires exchange of material between nanoparticles. By 
constructing phase diagrams for both the coherent and 
mosaic scenarios, we find a limited role for coherent nu- 
cleation and growth and predict that moderately ele- 
vated temperatures ought to suppress all two-phase co- 
existence, even in large particles. 
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Phase-field model 

We begin with the reaction-hmited phase-field model 
for ion-intercalation in single, anisotropic nanoparticles 
from our prior work (TTJ ISTJ |32] and extend it to in- 
clude coherency strain. The theory couples electrochem- 
ical surface reactions to bulk phase separation using a 
thermodynamically consistent generalization of Butler- 
Volmer kinetics. The reaction rate depends on the Cahn- 
Hilliard [33 (or Van der Waals [34 ) gradient energy, in- 
troduced to model the formation of phase boundaries. 

Individual nanoparticles are modeled as open thermo- 
dynamic systems in contact with an electrolyte mass 
reservoir at constant temperature, volume, and chemi- 
cal potential (grand canonical ensemble). For solids, PV 
work is generally very small and can be neglected at at- 
mospheric pressure. Thus we consider the free expansion 
of the particle at zero applied pressure. The inhomoge- 
neous grand free energy functional is [191 HO] • 

/I 1 
/(c) - /ic+ -A>:(Vc)^ + -Cijkieijeki dV 

(1) 

where c{x) is the mole fraction of lithium, /(c) is the 
homogeneous Helmholtz free energy density, /ic is a Leg- 
endre tranform that accounts for the chemical potential 
ji of lithium ions in the reservoir, k is the gradient en- 
ergy coefficient that introduces interfacial energy, and 
^Cijki^ij^ki is elastic strain energy. Cijki is the elastic 
stiffness tensor, and eijix) is the total strain field which 
may be decomposed into three parts: 



4c(f) 
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where Cij is a homogeneous strain that accounts for uni- 
form deformation of the particle, ^ {uij + Uj^i) is a local 
inhomogeneous correction to Cij resulting from composi- 
tional inhomogeneity, and e^j is the lattice misfit between 
FeP04 LiFeP04 stress-free strain). We assume 
that the misfit strain varies linearly with composition 
(Vegard's Law). 

The electrochemical reaction at the surface of the par- 
ticle is governed by a generalized Butler- Volmer equation 
for solids and concentrated solutions [35l [36] : 
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where the overpotential has a variational definition: 
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A<peq is the Nernst equilibrium potential and Acj) = — 
(p = ^ is the interfacial voltage, where (p and are the 
electrostatic potentials of ions and electrons, respectively. 
When a potential is applied, the system is displaced 
from equilibrium and lithium enters or leaves the system. 



It is a subtle but important point that we define over- 
potential relative to the Nernst voltage rather than the 
voltage plateau of the phase-separated system at zero 
current. The Nernst voltage is an equilibrium material 
property, but the voltage plateau is not. A fiat voltage 
plateau is commonly cited as a hallmark of phase separa- 
tion, but as we will demonstrate, coherency strain leads 
to upward-sloping plateaus. 

To model the experimentally relevant case of galvano- 
static discharge, the current ffow into the particle is con- 
strained by an integral over the active area at the surface 
of the particle: 



dc ^ . dc 



■e (5) 



where ^ is a Langevin noise term. Mechanical equilib- 
rium equations must additionally be solved for the dis- 
placement vector u{x): 
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Cijkieij{x) dV = (6) 



The second equation defines the average stress in the sys- 
tem to be zero, which is necessary for stress-free bound- 
aries and equivalent to minimizing S with respect to Cij. 

Small particles are expected to be limited by surface 
reactions, and in this case Eq. [5]and[6]constitute a depth- 
averaged system of equations that can be solved in 2D on 
the particle's active surface. The validity of the depth- 
averaged approximation is supported by anisotropic elas- 
tic considerations (Fig. [T]) and by phase-field simulation 
of diffusion in the lithium channels [37]. 



RESULTS 
Phase boundary orientation 

Since LiFeP04 is orthorhombic, it is necessary to 
consider its fully anisotropic elastic stiffness and the 
anisotropic lattice mismatch between phases when an- 
alyzing phase boundary morphology. Assuming that the 
elastic modulus of each phase is the same (homogeneous 
modulus assumption), Khachaturyan [38, 39 related the 
elastic energy of an arbitrarily anisotropic elastic inclu- 
sion to a function of direction: 



B{n) = Cijkie^ijeli - nia^jnji{n)af^nr, 



(7) 



n is the interface normal, and which is related to the 
elastic Green's tensor, is defined by its inverse tensor 
Q^j^ = CikijUkfii. Elastic energy is a function of orien- 
tation ft because a phase boundary produces zero strain 
in the normal direction. The direction no that minimizes 
Eq. [T] defines the habit plane, which is the elastically pre- 
ferred orientation of the phase boundary that minimizes 
strain energy. 

Equation [7] is plotted in Fig. [Taj using anisotropic 
elastic constants calculated via first-principles [4 and 



3 



[001] 



[001] 



Coherent 



Semicoherent 



[010] 




[101] 



[100] 



[010] 




[0 1] 



[100] 



/ 


I 1 

n i 


.5 .25\ 






(c) 



FIG. 1: Spherical plots of_B (n) , the elastic strain energy of a flat interface as a function of normal direction for (a)_ 
coherent interface and |(b)| a semicoherent interface which has lost coherency in the [001] direction. The red arrows 



indicate no, the direction of minimum energy, (c) A polar plot of B{n) in the in the a-c plane comparing coherent 

and semicoherent energies. Energy is in units of GPa. 



anisotropic lattice mismatch measured in [10] . The figure 
reveals that {010} and {001} interfaces are high-energy, 
which justifies a posteriori the depth- averaged approx- 
imation in our model [31 , based on fast diffusion and 
no phase separation in the {010} depth direction. The 
orientation of no obtained by numerical minimization is 
drawn in red. There are four minima which lie along the 
{101} family of crystal planes, where B{no) = .19 GPa. 
Fig. [2] compares the predicted {101} phase boundaries 
with experimental observations by Ramana et al. [12j . 
Laffont et al. also appear to have observed a {101} in- 
terface in Fig. 2b of [11 . 

On the other hand, {100} interfaces have been reported 
in some experiments [lOl H^, which according to Fig. 
[Tal should not be elastically preferred. The resolution 
to this apparent discrepancy may be the formation of 
dislocations, which lead to a loss of coherency in the [001] 
direction. Indeed, Chen et al. [W report observing cracks 
and dislocations running in the [001] direction of negative 
misfit strain. 

Stanton and Bazant [40 recognized the importance of 
negative misfit for LiFeP04 with isotropic elastic analy- 
sis, but here we consider the fully anisotropic case. Fig. 
lb plots B{n) for a semicoherent interface with 633 = 0, 
and Fig. [Tc] compares the coherent and semicoherent 
cases in cross section. The semicoherent habit plane lies 
along the {100} family of planes, and curiously B{no) re- 
mains nearly unchanged by to the loss of coherency. The 
orientation of the interface changes, but its elastic energy 
does not. 

The mechanism by which coherency is lost remains to 
be determined. It could be that the phases initially form 
coherently, but then lose coherency over time lose as dis- 
locations form. It is also possible that phase boundaries 
form semicoherently upon lithiation, aligned with pre- 
existing cracks or defects. The dynamics of phase sepa- 



ration for the latter scenario is presented in Fig. |6] 



Morphology 

Modulated structures resulting from coherency strain 
are often observed in experimental systems, and stripes 
are an equilibrium morphology that minimizes energy in 
finite size particles [39l [41] . Stripes form due to elastic 
relaxation at the surface of the particle, and align nor- 
mal to no- The characteristic wavelength A of the stripes 
balances the elastic energy of surface relaxation which 
scales with volume, and total interfacial energy which 
scales with particle size [39j |4T]. Evidence of this re- 
laxation is visible near the boundaries of the simulated 



particle in Fig. 2b 



The wavelength of periodicity is described by a scaling 
relation derived in the supporting material: 



A/ 



(8) 



A is the period of the striping, 7 is interfacial energy, Lc 
is the width of the particle in the [001] direction, and A/ 
is the difference in free energy density between the ho- 
mogeneous state and the coherent phase- separated state. 
A/ has a chemical contribution from the homogeneous 
free energy density f{c\ and an elastic contribution from 
coherency strain (Eq. [9|. We find A/ = 4.77 MJ/m^ us- 
ing the regular solution model and gradient energy that 
were fitted to experimental data in Methods. 

Phase-field simulation and experimental observation 
of stripes are compared in Fig. [2a| and [2c| Using 
Eq. [8| the striping in Fig |2d| can be used to obtain 
the FeP04/LiFeP04 interfacial energy. Applying Eq. 
|8] to the striped pattern in Fig. 2d, with Lc = 4/im 
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(a) 500x500 nm 



(b) 500x500 nm 




FIG. 2: Comparison of simulated and experimental 
microstructures in Li 5FeP04. See Fig. [6] for simulation 



dynamics. |(a)| Phase boundaries align alo ng { 101} 
planes to minimize elastic coherency strain, (b) Loss of 



coherency in the [0 01] d irection causes the stripes form 
along {100} planes. [(c)]HRTEM image of a {101} phase 
boundary. Repr inted from [12] with permission from 
Elsevier, [(d)] TEM image of {100} stripes [10 . 
Reproduced by permission of The Electrochemical 
Society. 



and A ~ 250 nm (measured away from the corner to 
mitigate the influence of particle geometry), we infer a 
phase boundary energy of 37mJ/m^. For the phase- 
field simulation in Fig. |2b[ Lc = 500 nm, A ^ 90 nm, 
and 7 = 39mJ/m^. The interfacial energy can also be 
calculated directly from the phase-field model [33] [42] . 
which yields 7 = 39mJ/m^, in agreement with the rule 
of thumb that coherent phase boundaries have interfa- 
cial energies less than 200mJ/m^ [43 . This confirms 
both the validity of the scaling relation and our choice of 
phase-field parameters (in particular which has until 
now been difficult to estimate). 

Fig. |6] shows the dynamics of phase separation for a 
homogeneous particle that is held at zero current. Both 
the coherent and semicoherent cases are shown. The ini- 
tial decomposition is followed by a period of coarsening, 
but coarsening stops when the stripes reach their char- 
acteristic wavelength, which scales with ^/L^. Thus the 
stripes are dependent on particle geometry and do not 
coarsen, as would be expected if they were related to the 
must unstable wavelength of spinodal decomposition ^12j . 
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FIG. 3: Comparison of calculated and measured 
LiFeP04 solubility limits as a function of temperature 
and particle size, (a) Size depende nce a t room 
temperature. Data points from [28 . |(b)| Size and 
temperature dependence. Data points from [9]. 



Critical particle size 

Phase-field methods have been used in studies of size- 
dependent solubility without coherency strain [36[ [44], 
and a minimum system size was found below which two- 
phase coexistence is prohibited. This minimum size is 
set by the diffuse width of the phase boundary. Here we 
find that the critical particle size criterion changes with 
the introduction of coherency strain, and is a result of 
the combined effects of coherency strain and interfacial 
width. 

Figure [3] compares phase-field calculations and mea- 
surements of the solubility limits as a function of particle 
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size and temperature. A regular solution was used for 
/(c), and the regular solution parameter ft and gradient 
energy were obtained with a least-squares regression of 
the phase-field model to the data points in Fig. [s] (see 
Methods). With just these two parameters we were able 
to simultaneously fit both size and temperature depen- 
dence of lithium solubility (four experimental data sets). 
This confirms that the Lia^FeP04 system may reasonably 
be described as a regular solution. 



The fitting in Fig. 3b offers new insight into the exper- 
iments themselves. In Ref. [9 , the differences in the mis- 
cibility gap were originally thought to be related to par- 
ticle size, but according to Fig. [Sj 100 nm and 42 nm par- 
ticles should not show significantly different solubilities. 
This assertion is confirmed by x-ray diffraction (XRD) 
analysis in the paper, which found strain in samples A 
and B, but not in C. Sample A and B were prepared dif- 
ferently than C, which may explain why phase boundaries 
did not form in the sample C particles (see below). Fig. 
3b shows that the shrinking miscibility gap in samples A 
and B is plausibly explained by coherency strain. 

The equilibrium phase boundary width was measured 
from simulation to be 12 nm, which is in good agreement 
with the 12-15 nm width measured by STEM/EELS [IT . 
In phase-field simulations of small particles near the crit- 
ical size, we observed that phase separation always oc- 
curs as a sandwich (see Fig. [3a| inset), sometimes with 
the lithiated phase in the middle, and other times with 
the delithiated phase in the middle. Presumably this is 
a result of elastic interaction between the phases. Both 
cases require the formation of two interfaces, explaining 
why the critical particle size of 22 nm is roughly twice the 
interfacial thickness. 



Phase diagram 

Fig. |4] shows a phase diagram that was calculated us- 
ing the fitted regular solution model (See Methods). A 
mosaic phase diagram was calculated using /(c), and a 
coherent phase diagram was calculated by adding elas- 
tic energy ( Eq. [9| to /(c). The eutectoid reaction 
[H [7] involving a disordered phase at higher tempera- 
tures has been neglected. The phase diagram reveals 
that coherency strain stabilizes the solid solution at tem- 
peratures above 150°C, well below the disordering tem- 
perature. 

The illustrations in Fig. [4] depict a completely delithi- 
ated cathode nanoparticle that has been discharged to 
the corresponding points in the phase diagram. At point 
A the particles are inside the mosaic miscibility gap, but 
do not transform since there is no phase transformation 
pathway. The microstructure is thus metastable with 
respect to mosaic decomposition. At point B, the parti- 
cles cross the coherent spinodal and coherent nucleation 
inside particles becomes possible. Phase transformation 
will proceed slowly in this region since nucleation is an 
activated process, and for fast discharge this region will 



be bypassed. 

By point C the particles have crossed the mosaic spin- 
odal, and spontaneously form a mosaic if they are able to 
exchange ions through the electrolyte (light blue). Cur- 
rent plays as important role in the onset of this transi- 
tion, since the position of the spinodal moves inward and 
eventually disappears with increasing current [17 . 

Point D illustrates the case where discharge occurs 
rapidly to a point inside the coherent spinodal without 
time for exchange of ions between particles. The red 
background color indicates that the particles are not in 
electrolyte, and therefore cannot exchange ions. The par- 
ticles relax to a phase-separated state in this case only. 
Chemical delithiation, which has produced observations 
of phase boundaries, corresponds to case D. 

The difference in solubilities predicted in Fig. |4] can 
be used as a guide to the interpretation of experimental 
data. If a mosaic forms, the existence of a second phase 
first becomes possible inside the mosaic solubility limits, 
and particles will have compositions of either = .01 
or X(3 = .99. XRD measurements of a system of fully 
intercalated and fully deintercalated individual particles 
confirms the appearance of a second phase by x = .04 at 
room temperature [21^, and equilibrium measurements 
find very little room temperature solubility [H [25] . 

However if coherent phase boundaries form inside par- 
ticles, the onset of phase separation in Fig. [4] will only 
occur for .09 < x < .91. Frequent reports of extended 
regions of solid solution [S^JTO^, '26'-'28^ support this 
prediction, and several authors [23^ 26 have attributed 
their observations to retained strain. Moreover, the XRD 
measurements of Chen et al. [lOj on striped particles 
did not detect phase-separation until at least 10% of the 
lithium had been extracted. Badi et al. ^\ also recently 
measured the composition of coexisting Lic^FeP04 and 
Li/3FeP04 phases, and found a ^ .1 and ^ ^ .9. Both 
of these observations agree precisely with the predicted 
coherent miscibility gap in Fig. [4] at room temperature. 



Phase separation dynamics 

Spinodal decomposition in systems with coherency 
strain was studied by Cahn [19, 20 , who found that for 
homogeneous systems at equilibrium, strain energy can 
be approximated for small fiuctuations as: 



T^CijkieijCki - ^B{no){c- x) 



(9) 



Since strain energy is a function of the mean composi- 
tion of the system, coherency strain invalidates the com- 
mon tangent construction and leads to an upward-sloping 
voltage plateau [24J. 

A linear stability analysis of the evolution equations 
(Eq. [5| was performed in [ITj, and the amplification 
factor was found to be: 
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FIG. 4: Phase diagram comparing coherent and mosaic phase separation. Sohd hnes are the hmits of miscibihty, 
and dashed hnes indicate the spinodal. Hashing denotes the region where coherent nucleation is feasible. The 
ihustrations depict a fuhy dehthiated cathode that has been discharged to points A,B,C and D in the phase 
diagram. A, B, and C are bathed in electroyte (hght blue), but D is not. 



Bar notation indicates evaluation of functions at the ho- 
mogeneous state c = X. The coherent diffusion potential 
is found by inserting Eq. |9] into Eq. [l] and taking the 
variational derivative: 

/i = /(c) - kV^c + 5(no)(c - x) (11) 

Linear stability with and without coherency strain is 
presented in Fig. [5a| and shows a transition from phase- 
separating to homogeneous filling as current increases. 
Coherency strain promotes the stability of the solid so- 
lution by reducing the critical current and shrinking the 
spinodal. The neutral stability curves in the figure repre- 
sent the boundary between stable and unstable dynamics 
(5 = 0), and the 5 = 1 curve is where the amplification 
factor is large enough to produce phase separation on the 
order of the discharge time. This curve is an indicator of 
when complete phase separation is observable, and has a 
maximum at ^ = .047. In between the coherent curve 
and the 5 = 1 curve is a region of quasi- solid- solution 
where there are unstable modes, but not enough time for 
complete phase separation [? ]. 

The transition from fully phase-separating to quasi- 
solid-solution is captured in the simulated microstruc- 



tures of Fig. |5]at x = 0.6 for filling at different currents. 
By = .05 (slightly above the maximum of the 5 = 1 
curve), phase separation is just barely visible. Thus we 
conclude that due to coherency strain, phase separation 
is suppressed when the applied current exceeds only a 
few percent of the exchange current. 



Voltage curves during discharge were calculated at dif- 
ferent currents and are presented in Fig. [Sj The most 
striking difference compared to the incoherent case [17\ 
is the upward-sloping voltage plateau when phase sepa- 
ration occurs, which was predicted at equilibrium by Van 
der Ven et al. [24 . The first ions to enter the particle do 
extra mechanical work straining the surrounding lattice, 
and this work is recovered by the last ions, which en- 
ter lattice sites that have already been partially strained. 
When phase separation is suppressed at higher currents 
however, there are no phase boundaries and hence co- 
herency strain does not play any role. 
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FIG. 5: Analysis and simulation of galvanostatic 
discharge, (a) Linear stability boundary as a function of 
Constant-current discharge curves for a 



current. 

single nanoparticle with coherent interfaces. |(c)|(e)| 
Phase boundary morphology showing formation of a 
quasi-solid-solution as discharge rate is increased 
(Li6FeP04, 100x100 nm). 



DISCUSSION 

In this paper we have presented a thermodynamically 
consistent phase-field model for nanoparticulate interca- 
lation materials and focused on the significance of co- 
herency strain in the Li^FeP04 two-phase system. With 



just two free parameters (the regular solution parameter 
Q and the gradient energy tv)^ our model simultaneously 
explains the observed phase boundary orientations, stripe 
morphologies, the measured phase boundary width, in- 
terfacial energy, size- and temperature-dependent solu- 
bilities, and reports of extended solid solution. Elastic 
analysis reveals that {101} is the preferred phase bound- 
ary, that negative misfit strain along the [001] axis ex- 
plains the observation of {100} phase boundaries, and 
that elastic relaxation at the surface of the particles is 
the origin of stripes. Analysis and simulation of galvano- 
static discharge shows that coherency strain significantly 
suppresses phase separation during discharge and leads 
to upward-sloping voltage curves for single particles be- 
low the critical current. 

Although phase boundaries have been observed exper- 
imentally, we conclude that most electrochemical data 
is inconsistent with phase boundaries forming inside 
LiFeP04 nanoparticles. The method of sample prepa- 
ration appears to strongly infiuence the feasibility of 
phase boundary formation. We suspect that chemi- 
cal delithiation may be partly responsible for the ob- 
servation of phase boundaries by preventing the ex- 
change of lithium between particles. Chemical delithia- 
tion [8l|l0l[l2l|26l|28], a high degree of Li antisite defects 
[27 , particle size, and the synthesis technique itself |9|[23] 
all appear to infiuence the formation of phase boundaries. 

According to Fig. [4j temperature may play an impor- 
tant role in improving battery performance. Most bat- 
tery research focuses on room temperature operation, but 
room temperature discharge passes through the hatched 
region of Fig. [4] where nucleation is possible and the 
coherent miscibility is crossed before the mosaic spin- 
odal. However for temperatures greater than 70° C, the 
situation reverses, and the mosaic instability occurs be- 
fore coherent nucleation. Therefore, moderately elevated 
temperatures could be useful for stabilizing the solid solu- 
tion. Homogeneous particles have advantages for battery 
performance since they have larger active area for inser- 
tion and do not waste energy forming phase boundaries. 
Homogeneous particles also avoid the internal stresses 
caused by phase boundaries, which benefits cycle life. 

By deriving a simple scaling relation, we have been 
able to estimate the phase boundary energy by measur- 
ing the stripe wavelength in an experimental micrograph. 
Our mathematical formula predicts 7 = 39 mJ in precise 
agreement with the calculated value from the phase-field 
model. The fact that 7 is small may have significant 
consequences for the role of nucleation at small currents. 
Cahn and Hilliard showed that nucleation and growth 
competes with spinodal decomposition near the limit of 
metastability [20^, ^45^ . In the case of Fig. [Ij the energy 
barrier for nucleation in the hashed region is important. 
Applying the approximation of Cahn and Hilliard, we 
estimate the room temperature energy of coherent ho- 
mogeneous nucleation to be 102 kT at the mosaic spin- 
odal. Although this energy is fairly large, heterogeneous 
nucleation at the particle surface is a more likely path- 
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FIG. 6: An initially homogeneous system (Li^FeP04, 



500x500 nm) decomposes into regions of high and low lithium 
concentration when held at zero current. |(a)| (d) The phase boundaries align along {101} planes to minimize elastic 
coherency strain. '(e)|(h) Loss of coherency in the [001] direction results produces {100} interfaces. 



way in small nanoparticles with a large surface to volume 
ratio. The heterogeneous nucleation barrier is likely to 
be only a fraction of the homogeneous barrier, placing 
heterogeneous nucleation well within the realm of kinetic 
relevance. Energy barriers must generally be less than 
78 kT for observable rates of nucleation [43]. Therefore 
our future work will focus on accurately calculating the 
critical nucleus and heterogeneous nucleation barrier en- 
ergy in order to understand the role of nucleation. 



METHODS 



Numerical methods 



Eq. [5] and [6] were solved in 2D as a coupled sys- 
tem of differential-algebraic equations using finite differ- 
ence methods on a square grid. The Matlab function 
ode 15s was used for time integration with the integral 
constraint implemented using a singular mass matrix. 
Zero-flux boundaries were applied to c{x) , and zero pres- 
sure boundaries were applied to u. At each timestep it 
was necessary to solve for c(x), Ui{x)^ '^ji^)^ ^4^-> ^i^ ^2, ^3, 
and 66. 



Phase-field parameters 

The inputs to the phase-fleld model that must be se- 
lected for the LiFeP04 system are the homogeneous free 
energy /(c), the gradient energy coefficient the elas- 
tic stiffness Cijku and the lattice mismatch e^^-. We use 
GGA+U first-principles calculations of Cijki for FeP04 
[4 , and lattice mismatch measured in [10]. For /(c) we 
assume regular solution model: 

VLif{c) = nc{l - c) + kT{c ln(c) + (1 - c) ln(l - c)) (12) 

where Vli is the volume per lithium atom in the solid. 
The regular solution parameter Q and the gradient en- 
ergy Hi were obtained with a least-squares regression of 
the phase-field model to experimental measurements of 
the miscibility gap [9, 28 , illustrated in Fig. [sj Phase- 
field simulations were performed by allowing a square 
particle at x = .5 to relax to equilibrium at zero cur- 
rent. The miscibilities were then found by taking the 
minimum and maximum compositions in the equilibrated 
microstructure. 

The best fit was achieved with Q = 115meV/Li 



(ll.lkJ/mol) and 



3.13 X lO^eV/m (5.02 



10 J/m). The root-mean-square error was 2.3%. Both 
values are close to those used by Tang et al. [16l |37] . 



Sample C in Fig. 3b had no measured strain and large 
particles that do not exhibit a size effect. Thus it is influ- 
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w 



w 



stripe scaling 




0.2 0.4 0.6 
X in LixFeP04 

(b) 



Li5FeP04 



FIG. 7: (a) Illustration of the stripe morphology in 
Elastic relaxation occurs in the shaded 



regions at the {001} surfaces of the particle, (b) The 
influence of coherent phase separation on free energy. 



enced by /(c) only. The fact that a particularly good flt 
is achieved with sample C supports the use of a regular 
solution for Li^FePO/, . 



Here we adapt the arguments of Khachaturyan [39l [41] 
to derive an expression for the period of striping in fl- 
nite size Li 5FeP04 particles. Illustrated in Fig. [7a| the 
stripes form to balance elastic relaxation at thejOOl} 
surfaces of the particle (energy /volume) and total inter- 
facial energy (energy/area). The change in energy due 
to elastic relaxation (A/ = 4.77 MJ/m^) is illustrated in 
Fig. [Tbj and is found with a common tangent construc- 
tion applied to /(c) + ^B{no){c — x)^, as described in 
[191 EQI m] • Assuming that the width of relaxation w at 
the boundary is comparable to the size of the stripes, the 
change in energy due to relaxation and creation of phase 
boundaries is: 



(13) 



The equilibrium stripe size will minimize AE^ and so we 



solve for w when = 0: 







dw 

Solving for the stripe period X = 2w^ we obtain: 



(14) 



X = 2w 



A/ 



(15) 



The dependence reveals that striping is an effect of 
flnite size domains. As the domain size approaches infln- 
ity, the equilibrium state approaches two inflnite domains 
separated by one interface. 
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